Driving quantized vortices with quantum vacuum fluctuations 
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We propose to use a rotating corrugated material plate in order to stir, through the Casimir-Polder 
interaction, quantized vortices in an harmonically trapped Bose- Einstein condensate. The emergence 
of such vortices within the condensate cannot be explained with a computation of the Casimir-Polder 
potential based on the pairwise summation approach or on the proximity force approximation. It 
thus appears as a genuine signature of non-trivial geometry effects on the electromagnetic vacuum 
fluctuations, which fully exploits the superfluid nature of the sample. In order to discuss quantita- 
tively the generation of Casimir-driven vortices, we derive an exact non-perturbative theory of the 
Casimir-Polder potential felt by the atoms in front of the grating. Our numerical results for a Rb 
condensate close to a Si grating show that the resulting quantum vacuum torque is strong enough 
to provide a contactless transfer of angular momentum to the condensate and generate quantized 
vortices under realistic experimental conditions at separation distances around 3 /im. 

PACS numbers: 03.75.Lm,12.20.Ds,34.35.+a 
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Since the advent of Bose-Einstein condensation [T| , the 
specific features of ultra-cold atomic samples have turned 
them into a versatile tool useful well-beyond the field of 
laser cooling and in particular for precision experiments. 
The extremely narrow momentum distribution of such 
samples can indeed be used in atom interferometric sys- 
tems serving as inertial sensors or atomic clocks [2|. Re- 
cently, Bose-Einstein condensates (BECs) have attracted 
considerable attention as promisingprobes of the atom- 
surface Casimir-Polder interaction [3, [j] . Casimir forces 
can indeed give rise to noticeable effects in the condensate 
dynamics, affecting for instance the stability of a mag- 
netically trapped sample [5|, or provoking quantum re- 
flections when the BEC approaches a surface sufficiently 
slowly Q. To date, the more accurate measurements of 
the Casimir-Polder force performed with a BEC [7| rely 
on the tiny correction to the center-of-mass oscillation 
frequency induced by a nearby material plate [8] . 

The experiments realized or proposed so far involve 
only the translational motion of the ultra-cold atom sam- 
ple. We propose instead to probe Casimir forces through 
the contactless transfer of angular momentum from a 
rotating grating to a harmonically trapped condensate. 
Quantum vacuum or Casimir torques have been pre- 
dicted to arise between two gratings when their corruga- 
tions are not aligned [9j or between two flat but bircfrin- 
gent plates [lfl|. However they have not been measured 
experimentally so far. In this Letter, we take advan- 
tage of the cold atoms technology and propose a vacuum 
torque between a rotating grating and a BEC resulting 
in the excitation of quantized vortices inside the BEC. 
This approach brings several benefits. It takes advan- 
tage of one of the most prominent features of superfluid- 
ity, namely the quantization of velocity circulation [ll| . 
In contrast to the experiments reported in Refs.Q, in our 
proposal the Casimir-Polder potential does not induce a 



small correction onto the atomic motion, it rather plays 
the central part in the stirring process: in the absence 
of Casimir effect, the ultra-cold sample would acquire no 
angular momentum. The emergence of quantized vortices 
is thus a genuine signature of the atom-surface interac- 
tion mediated by the quantum vacuum field. A fortunate 
feature of BECs is their ability to be efficiently stirred 
by a slightly anisotropic rotating potential, enabling the 
production of one-vortex states through a weak time- 
dependent contribution |l2{ . This property turns conden- 
sates into highly sensitive probes of rotating potentials. 
The purpose of this Letter is to show that, for realis- 
tic experimental parameters, the nucleation of quantized 
vortices can be driven by a rotating Casimir-Polder po- 
tential. 

In order to enhance the Casimir-induced anisotropy of 
the potential felt by the condensate, it is necessary to 
position it at a distance of the order of the corrugation 
period, which invalidates [13W15J ] the proximity force ap- 
proximation (PEA) |16j | and the pairwise summation ap- 
proach (PWS) 17] usually employed to derive Casimir 
forces. Within PEA, the Casimir-Polder potential is ap- 
proximated by the result for a plane interface, taken at 
the local atom-surface distance. As discussed below, it is 
precisely the departure from the PEA and PWS approx- 
imations that provides the dominant contribution to the 
time-dependent potential driving the vortices. In fact, 
neither the PEA nor the PWS treatments yield a suffi- 
cient torque to generate vortices in the sample. In this 
sense, the proposed setup probes non-trivial geometry ef- 
fects through the nucleation of quantized vortices, which 
contrasts to previous experiments [jj measuring Casimir 
forces in a simpler planar configuration. The Casimir- 
Polder potential for a grating could also be probed via 
Bragg spectroscopy of the energy spectrum [18| of a BEC, 
or by measuring its center-of-mass (dipole) lateral oscilla- 
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FIG. 1: Proposed experimental setup. 



tion frequency [14j . We explore here how such geometry 
effects can excite modes of the atomic-field presenting a 
non-zero angular momentum. 

The Casimir- Polder potential for a corrugated sur- 
face has been computed beyond the PWS and PFA 
treatments up to first order in the corrugation ampli- 
tude [Tj, [I3. We develop here a more general non- 
perturbative method based on the scattering approach 
to dispersive interactions [19|, |2(|, suitable to address 
the larger amplitudes considered here (see also [21j for 
a scalar field model) . This formalism, similar to the the- 
ory developed in Ref.[22| for treating surface-surface in- 
teractions, is built on the analysis of the scattering of 
vacuum field fluctuations between the ground-state atom 
and the grating. By taking the dispersive electromag- 
netic responses of both atom and material medium into 
account, it covers the entire range of distances, from the 
unretarded short-distance van der Waals regime to the 
Casimir-Polder large-distance asymptotic limit. For the 
relatively large distances considered here, we are closer 
to the Casimir-Polder limit, but our numerical results 
are nevertheless exact and take dispersion fully into ac- 
count. A detailed derivation of the dispersive potential 
is presented elsewhere [231 ] - 

We consider the following setup, sketched in Fig|T] 
It consists in a time-independent harmonic optical trap 
with cylindrical symmetry around the vertical axis de- 
noted as O z , the center of which is situated close to and 
below a plate rotating at an angular velocity f2 around 
this axis. The trap center lies on the rotation axis. The 
plate is assumed to have a square corrugation profile of 
period d, and to be sufficiently large to be treated as in- 
finite. The trap is loaded with a BEC preferably realized 
with an all-optical method [24j j in order to avoid spurious 
magnetic effects such as the induction of local charges on 
the rotating plate. We assume the sample to be well 
described within the mean- field approximation |l|. We 
denote Uo(r,z) — ^muj 2 j_(x 2 + y 2 ) + ^mujQ Z (z — z ) 2 
the sum of the gravitational and trapping potential and 
Ucp(r, z, t) the time-dependent Casimir-Polder potential 



induced by the plate, where Zq is the location of the trap 
center along the z-axis after taking into account the grav- 
itational sag. The trap lateral position is at x = y = 0. 
The transformation of the Hamiltonian into the rotating 

frame can be written as H = H — QL Z , where H and 

H are the laboratory and rotating frame Hamiltonians 
respectively, while L z is the vertical angular momentum. 
The second term in the r.h.s. favors the emergence of 
quantum states with an upward vertical angular momen- 
tum, which implies the presence of vortices in the con- 
sidered superfluid sample. 

We assume that the size of the condensate is suffi- 
ciently small to perform a quadratic approximation of the 
Casimir-Polder potential in the transverse coordinates 
with an effective frequency wcp evaluated at the cen- 
tral sample vertical position z$ : Ucp{x, z) ~ C/cp(0, z) + 
^muiQp(zo)x 2 . We have used the translational invariance 
of our setup along the axis O y . These assumptions limit 
the transverse and longitudinal radius of the sample to 
respectively a quarter of the corrugation period d and a 
fraction of the distance to the plate: 



R± < c± -, 



Rz < C/l Z , 



(1) 



with c±_ ~ I and c/i "C 1 . Considering a sample 
of small transverse size allows us to take advantage of 
the considerable theoretical and experimental work per- 
formed on the rotation of harmonically trapped conden- 



sates through time-dependent quadratic potentials [11 



Besides, in this configuration the PFA treatment predicts 
a Casimir-Polder potential which is time-independent in 
the sample region and thus unable to transfer any an- 
gular momentum to the condensate (the PFA potential 
is piece-wise constant in the transverse coordinates with 
jumps following the plate discontinuities.). 

The macroscopic sample wave-function then satisfies 
the time-dependent Gross-Pitaevskii equation (GPE) in 
the rotating frame [l[ 
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where the effective quadratic potential U(r, t) becomes 
anisotropic under the influence of the Casimir-Polder in- 
teraction 



U(r,t) 
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/2. In Eq.@, g is a constant 



related to the s-wave atomic scattering length a s by g 
4irh 2 a s /m. The anisotropy e is expressed as 



^CP 



/(2w, 



J CP 



(3) 



The normal Casimir-Polder potential Ucp{0, z) is not rel- 
evant in the following discussion. It slightly shifts the 



sample vertical position, which is taken into account by 
a redefinition of zq. After a stage of adiabatic approach 
of the plate, the rotating frame potential U(r, t) becomes 
time- independent . 

In order to optimize the sensitivity of our setup, we 
maximize the anisotropy e for a given plate corruga- 
tion and a given distance zq between the sample cen- 
ter and the plate. Accordingly, we choose the minimum 
radial trap frequency providing a confinement satisfy- 
ing Eqs.(fT]). Using the relations fj, = ^w^ z R±,z and 

|U = ^ftZJ [l5A^a/(/i/mw) 1 / 2 ] with ZZ? = u Q ' ± uj ' z be- 
tween the Thomas- Fermi radii R± yZ , chemical potential 
fi and trap frequencies uj±_ z , one finds readily 



u) x(d, z Q , N) = 



wVIEhiNas) 1 / 2 1 



1/2 2 

mc, Cj_ 



z d 2 



(4) 



This equation suggests to use a relatively dilute conden- 
sate, for which a weak radial trapping can balance the re- 
pulsive interactions and achieve the desired confinement. 
The pancake-shaped samples considered is this Letter 
are thus not rigorously in the Thomas-Fermi regime, i.e. 
their transverse radius R± is comparable to the the har- 
monic trap scale R u = y/h/(muj). However, the trans- 
verse kinetic energy is much smaller than the interaction 
energy, so that the required trapping frequency is given 
by Eq.((3]) with a good approximation. 

Let us now calculate the anisotropies that can be at- 
tained through the Casimir interaction between the con- 
densate and the grating. For the low angular frequen- 
cies used here, of the order of the transverse trapping 
frequency and thus of a few dozens of hertz, dynami- 
cal Casimir and non-contact quantum friction effects are 
negligible J25[. For separation distances in the microme- 
ter range, the time it takes for light to travel between the 
grating and the atom is of the order of 10 -14 s and hence 
much shorter than the time scale associated to the grat- 
ing rotation. Thus, the atom interacts with the instan- 
taneous angular position of the plate (on the other hand 
retardation is very important as far as charge fluctuations 
in both atom and grating are concerned). Therefore in 
the rotating frame the potential is given by the static 
potential, which we compute by following the scattering 
approach [13, [23[ . The atom-surface potential is written 
in terms of the dynamic atomic polarizability a evaluated 
at imaginary frequencies it; and of the reflection opera- 
tor IZs describing non-specular diffraction by the grating 
15]. 

To calculate TZs for a periodic grating, we employ 
the Rayleigh basis for the fields propagating inside the 
homogenous regions corresponding to the bulk material 
medium (z > 0) and to the empty space below the grating 
(z < —a). The fields in the modulated region of thickness 
a are obtained by solving coupled differential equations. 
The reflection matrix elements {j, a\TZs(k x , k v ,£)\j',a') 
(with k x varying in the first Brillouin zone [—ir/d, ir/d] 
and a and j representing polarization and diffraction or- 
der respectively) are then obtained by matching the ex- 



pansions in the three regions across the boundaries at 
z = and z = —a. The potential at position (x, y, z) is 
written as 



n r 00 d£ f°° dk y r /d dk x 

e c 2 J 2vr J_ 00 2?r J_„ /d 2tt 



p2 

x J2 —"(*0 e 27ri( - j - j ~> x ' d e -<*i+*i')H (5) 

xJ2(3,<r\Ks(k x ,k z ,O\j',a')^(k',iO-e+0t',iO, 



with k, 



k 2 . The unit vec- 



., - ^i 2 /c 2 + (k x + 2irj/dy 

tors e^r provide the direction of the electric field prop- 
agating upwards/downwards with polarization a = 
TE or TM. 

To calculate wcp, we first differentiate Ucp(x, z) twice 
with respect to x and then evaluate the resulting expres- 
sion at x — (middle of the plateau) numerically. We 
consider a silicon rectangular (lamellar) grating with am- 
plitude a = 4/im, period d = 30 /xm and gap width d/2 
(see Fig. 1). The corresponding polarizability function 
a(i£) is provided by Ref . [26| . whereas the electric per- 
mittivity for intrinsic silicon at the imaginary frequency 
axis, required to compute the reflection matrix TZs m ©, 
is calculated from data at real frequencies obtained from 
Ref. 12711 w ith the help of a suitable Kramers-Kronig re- 
lation [23|. We choose a sample of TV = 100 87 Rb atoms, 
which is both sufficient to detect vortices with state-of- 
the art techniques and compatible with a weak radial 
confinement. 

In Fig. [51 we plot the anisotropy e as a function of 
the separation distance, with ljo± computed from (|4|) 



with c± = 1 and c 



7/ 



0.2. We sketch the values ob- 



tained with the non-perturbative scattering method and 
with the standard PWS approach. These results can be 
compared with the regime of vortex nucleation in con- 
densates with quadratic rotating potentials, which have 
been extensively studied theoretically [llj and experi- 
mentally [29J, [30| . These studies determine a well-defined 
region in the parameter space of anisotropy e and angular 
frequency Q leading to the nucleation of vortices when 
a specific route is followed. Here we propose to follow 
the procedure of adiabatic increase of the anisotropy [29[ 
by slowly decreasing the distance between the rotating 
plate and the trapped sample until vortices emerge. Since 
the angular frequency of the plate can be tuned at will, 
we focus on the required anisotropy e which limits the 
separation distance at which vortices can be nucleated. 
Quantized vortices have been obtained experimentally 
with anisotropies as low as e mm = 2.23% [30(. Even if 
these studies have been conducted in the Thomas-Fermi 
regime and are thus not quantitatively applicable to the 
dilute samples considered here, they still provide quali- 
tative bounds which strongly support the emergence of 
vortices in the present setup. Fig.[2]shows indeed that the 
anisotropy induced by the Casimir-Polder interaction, as 
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FIG. 2: (Color online) Anisotropy induced by the exact (full 
red line) and PWS (blue dashed line) Casimir- Polder poten- 
tial as a function of the distance between the plate and the 
sample. We consider a silicon grating with period d = 30 /im 
and amplitude a — 4(j,m and a pancake-like condensate of 
N = 100 Rb atoms harmonically confined with a transverse 
radius R± — d/4 and a longitudinal radius of R z — 0.2\zq\. 
The horizontal dashed line represents the minimal anisotropy 
for which vortices have been produced experimentally. 



computed with the scattering method, is above the ex- 
perimental minimum e mm (dashed horizontal line) for 
distances up to 3 /xm. The pairwise approach predicts an 
anisotropy epws "well below the exact theoretical value 
e scat an d below the experimental threshold for vortex 
nucleation e mm . At a distance \zq\ = 3 /xm, one finds 
indeed epwsAscat ~ e pws/e m i n ^ 5 %. Thus, the emer- 
gence of vortices cannot be explained by the naive picture 
of a rotating additive atom-surface potential. 



To conclude, we have presented evidence that the elec- 
tromagnetic quantum vacuum fluctuations in the neigh- 
borhood of a rotating grating can transfer angular mo- 
mentum and induce quantum vortices in a condensate 
placed a few microns from the plate. In the setup we 
have analyzed, PFA would predict a fiat Casimir-Poldcr 
potential in the region of the condensate, whereas the 
PWS treatment predicts a distortion which is far too 
small to nucleate vortices in the sample. The appear- 
ance of vortices thus results from a non-trivial geometry 
effect on vacuum fluctuations, which is brought into evi- 
dence thanks to the supernuid nature of the sample. For 
the experimental parameters considered in our numerical 
example, significant anisotropics - higher than those re- 
quired to obtain vortex nucleation in the Thomas-Fermi 
regime - can be achieved for distances up to 3 /im. One 
can also explore the quantum state engineering of other 
BEC modes by tayloring quantum vacuum fluctuations 
with different non-trivial plate geometries. Last, we un- 
derline that we have proposed the quantum analogue of 
a rotating bucket which does not touch the sample: this 
setup illustrates, through the nucleation of quantized vor- 
tices, how quantum fluctuations may mimic contact in- 
teractions at a distance and provide contactless transfer 
of angular momentum. 
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